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ABSTRACT 



The imaging spectrometer SPI on board ESA's INTE- 
GRAL observatory provides us with an unprecedented 
view of positron annihilation in our Galaxy. The first 
sky maps in the 511 keV annihilation line and in the 
positronium continuum from SPI showed a puzzling con- 
centration of annihilation radiation in the Galactic bulge 
region. By now, more than twice as many INTEGRAL 
observations are available, offering new clues to the ori- 
gin of Galactic positrons. We present the current status 
of our analyses of this augmented data set. We now de- 
tect significant emission from outside the Galactic bulge 
region. The 511 keV line is clearly detected from the 
Galactic disk; in addition, there is a tantalizing hint at 
possible halo-like emission. The available data do not yet 
permit to discern whether the emission around the bulge 
region originates from a halo-like component or from a 
disk component that is very extended in latitude. 

Key words: Galactic 511 keV line emission; positron an- 
nihilation. 



1. INTRODUCTION 



The detection of 5 1 1 keV positron annihilation line emis- 
sion from the central region of our Galaxy was one of 
the early and important successes of gamma-ray astron- 
omy 11171 12211 . Although positron annihilation gives rise 
to the strongest gamma-ray line signal from our Galaxy, 
three decades of dedicated observational and theoretical 
effort have fallen short of unveiling the origin of Galac- 
tic positrons. A large variety of potential astrophysical 
and exotic positron sources have been proposed, includ- 
ing: interactions of cosmic rays in the interstellar medium 
B32I1 ; pulsars 13911 and millisecond pulsars Il44ll : the de- 
cay of radioactive nuclei expelled by supernovae |0], no- 
vae 0], Wolf-Rayet stars j|29ll . or hypernovae/gamma- 



ray bursts fH]; c omp act objects harbouring black holes 
or neutron stars (33ll. low-mass X-ray binaries lf30ll . and 
micro-quasars 1121 : gamma-ray bursts 12311 : pair produc- 
tion by gamma rays from so-called Small Mass Black 
Holes with X-ray photons from the supermassive black 
hole at the center of our Galaxy, Sgr A* l4lll ; interactions 
of cosmic rays accelerated during episodes of mass ac- 
cretion onto the supermassive black hole Sgr A*J5lj421; 
and decay or annihilation of (light) dark matter 11351 |2I1 or 
other exotic particles and processes [see e.g. references 
in Ql. 

Investigations of the sky distribution of the annihilation 
radiation promise to provide clues for the identification of 
the source(s) of positrons in our Galaxy. Positrons may 
travel from their birth places before annihilating, how- 
ever, recent theoretical work suggests that positron diffu- 
sion is limited 11131 flOil . Thus the Galactic distribution of 
annihilation radiation can be expected to resemble that of 
the positron sources. 

First maps of the annihilation radiation, limited to the in- 
ner regions of our Galaxy, were obtained using the OSSE 
instrument on board the Compton Gamma-Ray Observa- 
tory in the 511 keV line and in positronium continuum 
emission [e.g. ED, H |H |H H| H3]. Some of these 
analyses also took advantage of data obtained with the 
Gamma-Ray Spectrometer on board the Solar Maximum 
Mission (SMM) and the TGRS spectrometer on board 
the Wind satellite. Furthermore, the OSSE instrument al- 
lowed Kinzer, Purcell, & Kurfess lfl8ll and Kinzer et al. 
lfl9ll to study the one-dimensional distribution in longi- 
tude and in latitude of diffuse emission, including annihi- 
lation radiation, from the inner Galaxy. 

With the commissioning of the imaging spectrometer SPI 
on board ESA's INTEGRAL observatory, high spectral 
resolution mapping with improved angular resolution has 
become feasible JlllSIIalMllTl]. Both the 511 keV 
line emission and the positronium continuum emission 
were found to be brightest in a region of a few degrees 
around the Galactic center; emission from the Galac- 



tic disk is much fainter, implying that positron annihi- 
lation is concentrated in the central regions of our Galaxy 
llHil. In particular, the high bulge-to-disk flux 
ratio of 1-3 (corresponding to a luminosity ratio of 3- 
9) measured in the annihilation line imposes severe con- 
straints on potential positron sources. The most promis- 
ing positron sources for the Galactic bulge were found 
to be members of the old stellar population, specifically 
Type la supernovae (SN la) and/or low-mass X-ray bina- 
ries (LMXBs). Annihilation or decay of light dark matter 
was another possibility. At least part of the faint emission 
from the Galactic disk must be due to the (3 + -decay of the 
radioisotopes 26 Al and 44 Ti which are produced by mas- 
sive members of the young stellar population. However, 
none of these positron sources is without caveats [see e.g. 

Additional insight into the origin, propagation, and anni- 
hilation of positrons can be obtained from high-resolution 
spectroscopy. The detailed shape of the annihilation line 
and the ratio of the fluxes in the line and the positronium 
continuum depend on the physical conditions of the inter- 
stellar medium in which positrons annihilate [see e.g. fTlfl . 
Using observations of the Galactic bulge, Jean et al. IU5I1 
showed that positrons annihilate in the warm phase of 
the interstellar medium. Recently, Beacom & Yiiksel [1] 
pointed out that spectroscopy of Galactic emission above 
511 keV can be used to set an upper limit on the initial 
energy at which positrons are injected into the interstellar 
medium [see also [37tl . 

In this contribution, we present the current status of our 
studies of the large-scale distribution of Galactic 511 ke V 
annihilation line radiation using more than 2 years of ob- 
servations with the spectrometer SPI. An update on SPI 
spectroscopy of positron annihilation radiation is given in 
a companion paper by Jean et al. Biol . After describing 
our analysis methods, we first present an updated sky map 
in the 511 keV line. We then summarize our current re- 
sults regarding the distribution of the bright 511 keV line 
emission from the Galactic center region and regarding 
the existence and distribution of more extended emission 
from the Galactic disk and beyond. 



2. INSTRUMENT AND DATA ANALYSIS 



The SPI imaging spectrometer consists of an array of 19 
actively cooled high resolution Ge detectors, which cover 
an energy range of 20-8000 keV with an energy resolu- 
tion of about 2.1 keV FHWM at 511 keV. SPI employs 
an active anti-coincidence shield made of bismuth ger- 
manate (BGO), which also acts as a collimator. In addi- 
tion to its spectroscopic capabilities, SPI can image the 
sky with moderate spatial resolution of about 3° FWHM 
using a tungsten coded aperture mask. The fully coded 
field-of-view of the instrument is about 16°. A detailed 
description of the instrument was given by Vedrenne et 
al. El. 

The results presented here are based on a data set that 




Figure 1. Top panel: the exposure to the sky of the data 
set used in this analysis. Contour levels are at 1, 2, 3, 
and 4 X 10 7 cm 2 s. Bottom panel: the exposure to the 
sky after the first year of the mission, depicted with the 
same color coding than above. Contour levels are at 1 
and 2 x 10 7 cm 2 s. 



comprises more than twice as many observations as that 
we used in our studies after the first year of the mission 
llH [H IH. Specifically, we have used the April 20, 
2006 public INTEGRAL data release (i.e. three-day or- 
bital revolutions 16-269, 277, 278, 283-285). These pub- 
lic data were supplemented by instrument team observa- 
tions of the Galactic center region and of the Galactic 
plane and by recent private observations of SNR 1006, 
Cen X-4, and Sco X-l up to revolution 423. The obser- 
vations were taken during the epoch November 23, 2002 
through April 1, 2006. During this time a few very strong 
(notably during October 2003 and September 2005), and 
several smaller, solar flares occurred. These resulted in 
strong transient enhancements of the instrumental back- 
ground at 5 1 1 keV. These exceptional periods were re- 
moved, resulting in a cleaned data set that consists of 
18101 pointings with a combined live time of 3.6 x 10 7 s. 

The resulting exposure to the sky is depicted in the top 
panel of Fig. [T] for comparison, the exposure from the 
observations of the first year of the mission is shown in 
the bottom panel. To a first approximation, the exposure 
of the current data set is about twice that achieved af- 
ter the first year. Of particular importance for the anal- 
yses presented here is the fact that the exposure along 
the Galactic plane not only deepened, but also became 
more homogeneous within about \l \ < 90°. However, un- 
fortunately the substantial inhomogeneity outside that re- 
gion has not been alleviated. The exposure is still deepest 
around the central region of our Galaxy, and most of the 
remaining, lower, exposure away from that region is still 
concentrated in a rather narrow band of ±10° along the 
inner part of the Galactic disk. Investigations of the large- 



scale distribution of 511 keV line annihilation therefore 
still face the dilemma that regions of putative faint and 
extended emission remain poorly exposed. Specifically, 
the current inhomogeneity of exposure severely limits 
our ability to distinguish between models for low surface 
brightness emission outside the bulge region. This is par- 
ticularly pertinent to investigating a potential halo distri- 
bution around the Galactic bulge, and to determining the 
latitude distribution of the disk emission. 

In our analyses, we followed the procedures for mod- 
elling the instrumental background that we established 
earlier [for more details, see e.g. I20I1 . The count rate 
in the detectors is the sum of time variable instrumental 
background and signal from celestial sources. The latter 
is variable in time because even if the celestial source is 
intrinsically stable SPI's exposure to it varies in time as 
the instrument performs a series of observations. Source 
components have therefore to be distinguished from the 
dominant background components by taking advantage 
of their differing time variations. Hence our approach 
for extracting the source components consists of fitting 
time series of detector count rates by a linear combination 
of so-called background templates (time series of back- 
ground components, explained below) and mask modu- 
lated signals. 

We performed our analyses in a 5 keV wide energy in- 
terval centered at 5 1 1 keVfl In this energy range the in- 
strumental background consists of two components: the 
instrumental 511 keV back grou nd line, and an underly- 
ing continuum background 

TUBUS]. 

The variations 

in time of these two background components differ, so 
they are modelled independently ll2"oll . The time varia- 
tion of the 511 keV background line is modelled using 
a linear combination of three templates, the rate of satu- 
rating (> 8 MeV) events in the Ge detectors (GeDsat), 
a constant, and a linear increase in time. In this con- 
text, a template represents a specific time series of back- 
ground count rates in a given detector. The role of the first 
two templates is to account for prompt and short-lived 
background components, the third template accounts for 
long-lived background components. The time variation 
of the continuum background underlying the instrumen- 
tal 511 keV line is modelled by assuming it is identical 
(per unit energy) to that observed in an adjacent energy 
interval, for which we chose 523-545 keV. To reduce sta- 
tistical uncertainties, the continuum time variation was 
smoothed [see 12011 . 

When applying this 4-component background model, 
the normalizations of the three line background compo- 
nents are adjusted independently for each detector to ac- 
count for detector-to-detector variations of the instrumen- 
tal background. In addition, the normalizations of the 



'When varying the width of the analysis energy interval, the bulg e 
flux varies as expected from the spectral results by Jean et al. flal . 
Based on their two-Gaussian model for the 511 keV line spectrum, the 
flux determined in our 5 keV analysis interval corresponds to 82% of 
the total line flux. The variation of the disk flux is consistent with that 
of the bulge flux; detailed spectroscopy results from the current data set 
are reported in Jean et al. Ilql . 



GeDsat component are adjusted for each orbital revolu- 
tion, as in model ORBIT-DETE described in Knodlseder 
et al. i20ll . An additional complication arises from the 
fact that during the epoch covered by our data set, two SPI 
detectors failed (detector 2 on Dec. 7, 2003, or IJD 1435, 
and detector 17 on Jul. 17, 2004, or IJD 1659). Each time 
a detector fails, the detector-to-detector ratio of the in- 
strumental background changes significantly, especially 
in the 511 ke V background line and for detectors next 
to the failed detector. We therefore split the data set 
into three periods, each corresponding to a stable detec- 
tor configuration, and simultaneously apply three sepa- 
rate ORBIT-DETE models. 



3. RESULTS 



3.1. Imaging 

To obtain a model independent sky map of the 511 keV 
positron annihilation line radiation, we employed the 
Multi-Resolution Expectation Maximization (MREM) 
algorithm described in Knodlseder et al. 112111 . This al- 
gorithm is an extension of the implementation of the 
Richardson-Lucy (RL) algorithm that we applied in ear- 
lier analyses [e.g. KOIl . 

The resulting sky map is depicted in Fig. [2] As in ear- 
lier SPI analyses of the positron annihilation line ll20ll 
and the positronium continuum M47I1 . the only prominent 
511 keV line signal is that seen from the central region 
of our Galaxy. Any emission from other sky regions is 
much fainter. We are investigating the effect of the noise 
filtering threshold on MREM maps, as well as the effect 
of smoothing scales on RL maps, in an effort to deter- 
mine the best approach for uncovering faint and extended 
annihilation emission that we detect by model fitting (see 
Sec.0. 

As with the first year observations, we find again that the 
centroid of the emission is close to the Galactic center. 
The dominant emission from the central Galaxy appears 
to be more concentrated than seen in the first year ob- 
servations. This conclusion is confirmed by model fit- 
ting (see Sec. 13. 2\ . The total flux in the map is about 
1.0 x 10" 3 phcm- 2 s" 1 . 



3.2. Model Fitting 

A more quantitative approach for studying the Galac- 
tic distribution of the observed extended line emission 
is model fitting, which we performed using a maximum 
likelihood multi-component fitting algorithm described 
in Knodlseder et al. 112011 . In the present work, our focus 
is on investigating whether there is 5 1 1 keV line emission 
outside the central region of the Galaxy. In particular, we 
are trying to characterize the emission from the Galactic 
disk, and to assess whether there is evidence for emission 



Figure 2. An MREM sky map of the 511 keV positron annihilation line emission. The contours indicate intensity levels of 
10 -2 , 10~ 3 , and 10~ 4 /?/j era' 2 s^ 1 sr" 1 . Details are given in the text. 



from a halo component. In the first year data, disk emis- 
sion was already marginally detected, while a stellar halo 
component (comprising emission peaking at the Galac- 
tic center and fainter emission extending far beyond the 
bulge region) could not yet be discerned from pure bulge 
models. 

In order to investigate the possible existence of faint and 
extended emission from outside the central region of the 
Galaxy we fitted simple and flexible models for bulge 
and halo emissions in the absence or presence of simple 
disk models. The bulge and halos emissions were repre- 
sented by nested spherical shells of homogeneous emis- 
sivity, centered at the Galactic center. The disk emission 
was described either by the young (0-0. 15 Gyr) or the old 
(7-10 Gyr) stellar disk models as derived by Robin et al. 
I134ri - The fit results are summarized in Table Q] 

It is evident from Table[T]that there is significant 511 keV 
line emission from outside the bulge region of the Galaxy. 
The maximum likelihood ratio A for the "2 nested shells" 
bulge model is only 1273.8, much lower than for any 
other model tested. Adding a halo component (i.e. shells 
extending farther from the Galactic center) and/or a disk 
model significantly improves the fits. The disk mod- 
els are favoured over the halo models since either disk 
model improves the fit more than additional shell/halo 
components. In any combination of halo and disk com- 
ponents, the disk component is always significantly de- 
tected. When combined with a disk, the 1 .5-5 kpc shell 
is still marginally detected, which provides a tantaliz- 
ing hint at possible halo-like emission; the 5-8 kpc shell 
is then not required. The fluxes from the two inner- 
most shells, describing the bulge region, are remarkably 
robust and independent of the presence of other model 
components, reflecting the brightness of the bulge re- 
gion of our Galaxy in annihilation radiation. For the 
two-shell "bulge only" model, we obtain a total flux of 
(0.96 ± 0.06) x 10~ 3 ph cm" 2 s _1 , in excellent agree- 
ment with Knodlseder et al. j20Tl . For all other mod- 



els in Table [T] we obtain total fluxes in the range (1.7— 
3.1) xl0~ 3 ph cm~ 2 s _1 . Considering only models in- 
cluding disk components, but excluding models with four 
shells (since the large flux attributed to the 5-8 kpc shell 
is very uncertain), we obtain bulge+halo fluxes in the 
range (0.8-1.3) xl0 _3 ph cm -2 s _1 and disk fluxes in 
the range (0.9-1.8) xl0~ 3 ph cm~ 2 s _1 . The bulge-to- 
disk (B/D) flux ratio is found to be between 0.4 and 1 .4. 
These values are lower than the range 1-3 determined by 
Knodlseder et al. |20] using the first year of SPI obser- 
vations. Although mostly consistent within statistical er- 
rors, we tend to find now lower bulge fluxes and higher 
disk fluxes than with the first year data. Our values for 
the B/D flux ratio lie within the wide range of 0.2-3.3 
obtained by Milne et al. fH from OSSE/SMM/TGRS 
observations. 

We then investigated other characterizations of the bulge 
and disk emission. In a first step, and to compare with 
previous work, we modelled the bulge emission using an 
ellipsoidal distribution with a Gaussian radial profile in 
longitude and latitude, with full-widths at half maximum 
in longitude and latitude Ti and Tf,. This bulge model 
was combined with the young stellar disk description by 
Robin et al. |H; its three parameters are the disk scale 
length hp + , the scale length of the central disk hole ho , 
and the axis ratio e [see definitions in Table 3 of 13411 . 
We find values of T t = 6.5°±Hl and r b = 5.1°±£jo 
for the bulge component. These values indicate a bulge 
extent that is slightly smaller than the FWHM of about 
8° inferred from the first year data 12011 . We note that 
Kinzer et al. [H found values of Ti = 6.3° ± 1.5° 
and r 6 = 4.9° ± 0.7° in an analysis of OSSE ob- 
servations, in very good agreement with our measure- 
ment. The parameter values for the disk component are 
h>R + ^4 kpc, Iir_ ~ 3 kpc, and e ~ 0.3. As will be dis- 
cussed below, none of these disk parameters is well con- 
strained by the current data (for comparison, the young 
stellar disk is defined by Hr, — 5 kpc, hp_ — 3 kpc, 
and e = 0.014). Our disk parameters correspond to 



Bulge/Halo Model Bulge/Halo Flux Old Stellar Disk Flux A 

[phcm _2 s _1 ] [phcm -2 s -1 ] 

Dwek et al. (1995) stellar bulge E3 (8.67 ±0.38) x 1CT 4 (1.62 ±0.29) x 10" 3 1281.1 

Dwek et al. (1995) stellar bulge G3 (8.48 ± 0.36) x 10~ 4 (1.60 ± 0.29) x 10~ 3 1295.4 

Freudenreich (1998) stellar bulge S (8.87 ± 0.39) x 10~ 4 (1.57 ± 0.30) x 10" 3 1262.7 

Freudenreich (1998) stellar bulge E (8.87 ± 0.39) x 10~ 4 (1.56 ± 0.30) x 10" 3 1264.1 

Robin et al. (2003) stellar halo (2.35 ± 0.10) x 10~ 3 (1.16 ± 0.30) x 10~ 3 1319.6 

Table 2. The results of fits to the 511 keV line data of the preferred triaxial or bar- shaped stellar bulge models E3 and G3 
by Dwek et al. H^j and S and E by Freudenreich [9], and of the stellar halo model by Robin et al. [34]. A is the maximum 
likelihood ratio of the fits. 



a FWHM in longitude of about 70° for the 511 keV 
line emission, which is about twice as large as found by 
Kinzer et al. lfl9ll with OSSE. The maximum likelihood 
ratio for this model is A = 1315.4, and the bulge and 
disk fluxes are (7.04 ± 0.32) x 10" 4 ph cm~ 2 s" 1 and 
(1.41 ± 0.17) x 10~ 3 ph cm~ 2 s — x , respectively. In this 
case, the B/D flux ratio is about 0.5. Considering the 
number of model parameters that have been optimized, 
this description of the data using an ellipsoidal bulge and 
and a very extended (in latitude) disk is not statistically 
superior to the two-shell bulge models combined with a 
stellar disk (compare Table [T]). 

Rather than changing the aspect ratio of the Gaussian, 
more improvement is possible by changing the profile. 
The bulge model can be improved by using e.g. two 
Gaussian distributions instead of one ellipsoidal distri- 
bution. The values for the FWHM of the two Gaus- 

—1—1 A ^ — i-Q 

sians are 2.1°_ 1 ' 5 o and 8. O^j/qo, respectively. The value 
of A rises to about 1325, indicating that the 511 keV 
emission may be peaked or cusped at the Galactic cen- 
ter. The fluxes in the narrow and wide Gaussian and 
in the disk model are about 1.5 x 10~ 4 ph cm s , 
7.4 x 10~ 4 ph cirr 2 s" 1 , and 1.3 x 10~ 3 ph cm" 2 s -1 , 
respectively. Now the B/D flux ratio is about 0.7. 

To further investigate the distribution of the dominant an- 
nihilation radiation from the central region of the Galaxy, 
we also fitted to the 511 keV data some astrophysical 
models, namely the preferred triaxial or bar-shaped stel- 
lar bulge models E3 and G3 by Dwek et al. [8] and S 
and E by Freudenreich [9], and the stellar halo model by 
Robin et al. [ 34], each combined with the old stellar disk. 
The results are summarized in Table [2] The stellar bar 
models E3, G3, S, and E provide rather poor descriptions 
of the emission from the bulge region, compared to the 
simple shell models summarized in Table [T] and com- 
pared to ellipsoidal or Gaussian descriptions. It therefore 
appears that the distribution of the 511 ke V line emission 
does not follow the stellar bulge. The stellar halo model 
of Robin et al. Q, however, provides a fit that is of simi- 
lar quality than that of the best analytical models. As will 
be discussed below (see SecHJ), it is too early to conclude 
that indeed halo emission around the Galactic center re- 
gion exists. The good agreement between data and stellar 
halo model may mainly be due to the fact that it models 
well the peaked emission at the Galactic center. The ex- 



istence of a sharp peak in the emission is already hinted 
at by the good fit provided by the two-Gaussian bulge 
model described above. For the stellar halo model (which 
includes bulge emission), the total 511 ke V line flux is 
about 2.35 x 10~ 4 ph cm~ 2 s _1 , and the (bulge+halo)- 
to-disk flux ratio is about 2. This is consistent with the 
earlier analysis of Knodlseder et al. B20I1 and also with 
estimates of halo fluxes by OSSE/SMM/TGRS HIH. 

Our final step to date in investigating whether there is 
511 keV line emission outside the central region of the 
Galaxy, and whether it can be uniquely described by 
rather simple models, was to vary the axis ratio param- 
eter in the Robin et al. 113411 young and old disk models. 
Combined with either the two-shell (see Table[TJ, with an 
ellipsoidal (T; = 6.5° and = 6.0°) description of the 
bul ge e mission, or with the stellar halo model by Robin et 
al. Il34ll . As an example, for the young disk the variation 
of the maximum likelihood ratio A with the axis ratio e for 
the three descriptions of the emission from the central re- 
gion of the Galaxy are depicted in Fig. [3] The results for 
the old stellar disk are similar. We note that the axis ratios 
for the Robin et al. JH young (0-0.15 Gyr) and old (7- 
10 Gyr) stellar disks are 0.014 and 0.0791, respectively. 
For both bulge models, the preferred axis ratio is about 
0.3, much larger than the young disk value, in particular. 
The uncertainty is, however, rather large. For the stel- 
lar halo, the preferred value of e is similar to those for the 
Robin et al. [34] stellar disks, and the uncertainty is much 
smaller. These results clearly demonstrate that there is 
faint, extended emission around the central region of the 
Galaxy, but with the current data it is impossible to decide 
whether this emission is due to a halo-like component, or 
a disk very extended in latitude. 

As mentioned above, we have begun to study the longi- 
tude distribution of the Galactic 511 keV line emission 
by determining the best-fit values not only for the axis ra- 
tio e, but also for the disk scale length Iir + and for the 
scale length of the central hole hu_ for the Robin et al. 
ll34ll young stellar disk parameterization. Their best fit 
parameters correspond to a FWHM in longitude of about 
70° for the 511 ke V line emission. We have also begun to 
study the longitude distribution by dividing the disk emis- 
sion into longitude bins. First results indicate that we do 
not yet detect 511 keV emission from disk regions that 
are more than 50° in longtiude from the Galactic center. 
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Figure 3. Left panel: the variation of the maximum likelihood ratio A when fitting the data with the Robin et al. [34 ] 
young disk plus an ellipsoidal model (FWHM Fi = 6.5° and Tf, = 6.0°) of the bulge emission. Center panel: as left 
panel, except that the bulge emission is modelled by a combination of 0-0.5 kpc and 0.5-1.5 kpc homogeneous shells. 
Right panel: as left and center panels, except that the emission from the central region of the Galaxy is described by the 
stellar halo model of Robin et al. [34]. The dashed lines indicate decreases of X by 1 and 4, corresponding to 1 and 2a 
confidence levels for 1 degree of freedom. 



Both results indicate that 511 keV line emission is bright- 
est from the inner Galaxy. More exposure, in particular 
in the outer Galaxy, will be required to better determine 
the longitude profile. 



4. SUMMARY AND DISCUSSION 



The analyses we have performed to date, taking advan- 
tage of more than 2 years of observations, clearly demon- 
strate that SPI detects 511 keV annihilation line emission 
from outside the central or bulge region of the Galaxy. 
In similar, earlier analyses using only the first year of 
observations, the disk was only marginally detected by 
SPI, while the emission from the central region of the 
Galaxy could statistically be equally well described by 
pure bulge or bulge and halo distributions. Before the 
launch of INTEGRAL, positron annihilation from the 
Galactic disk had already been detected using OSSE ob- 
servations, sometimes supplemented by SMM and TGRS 
data [e.g.|2J,|l9|l- However, constraining the spatial char- 
acteristics of the disk emission, in particular its latitude 
distribution, eluded OSSE/SMM/TGRS, as did a conclu- 
sive result regarding the existence of halo emission. 

While SPI has not yet solved these two major goals 
of positron astronomy, the current status of our studies 
with SPI does represent significant progress and promises 
much more. SPI now detects disk emission at a signifi- 
cance level of about 6a (see e.g. Tabled, thus confirming 
the earlier detections with the OSSE/SMM/TGRS instru- 
ments. In addition, with SPI we find tantalizing hints at 
possible halo-like emission. When representing the halo- 
like emission by extended shells as in Table Q] we detect 
the 1.5-5 kpc shell at the ~ 2.5a level even when a disk 
model is included in the fit. Another indication for po- 



tential halo emission is the good fit of the emission from 
the central region of the Galaxy provided by the stellar 
halo model of Robin et al. 13411 . which comprises emis- 
sion peaking at the Galactic center and fainter emission 
extending far beyond the bulge region. Finally, the ex- 
istence of faint, extended emission around the Galactic 
bulge region is demonstrated by the fact that when a halo 
component is not included in the model disks with very 
wide latitude extents are favoured. When describing the 
central emission with a stellar halo, the preferred latitude 
extent of the disk is much smaller (see Fig. [3}. 

Despite this tantalizing indications for halo-like emis- 
sion, some caution is in order because of the inhomoge- 
neous exposure to the sky (see Fig. [TJ. With the current 
exposure to the sky, we cannot exclude that most or all 
of the flux attributed to the 1.5-5 kpc shell model is ac- 
tually originating from the Galactic disk, particularly in 
case the disk is extended in latitude. Similarly, the stellar 
halo model by Robin et al. [34] may describe the emission 
from the central region of the Galaxy well mainly because 
the model is rather peaked or cusped at the Galactic cen- 
ter. The existence of such a peak or cusp is indicated by 
the good fit provided by the two-Gaussian (FWHM about 
2.1° and 8.0°) bulge model. Since both the emission and 
the exposure are maximal around the Galactic center, it 
is possible that the fits result is driven by a rather small 
sky area in the bulge. Lower surface brightness emis- 
sion around the bulge would fall in sky regions that have 
received relatively little exposure. They therefore carry 
comparatively little weight in the fitting. 

More observations around the Galactic bulge region will 
be needed to alleviate the existing systematic biases due 
to the current inhomogeneity of the exposure. With the 
existing data, we cannot yet discern whether the emission 
around the bulge region originates from a halo-like com- 



ponent or from a disk component that is very extended 
in latitude. In particular, the latitude distribution of the 
disk cannot be determined independently of assumptions 
about the distribution of emission in and around the bulge 
region. Thus the situation reached with the INTEGRAL 
data with more than 2 years of data already starts to im- 
prove on what was possible with the combination of data 
from OSSE/SMM/TGRS at the end of the 9 year CGRO 
mission S HI- It is hoped that INTEGRAL will con- 
tinue to operate for many years and make further ad- 
vances possible. 



The current data indicate that the disk is brightest in 
511 ke V line emission in the inner Galaxy. Preliminary 
results indicate a FWHM in longitude of about 70° for 
the 511 keV line emission, which is about twice as large 
as determined by Kinzer et al. [ 19] with OSSE. More ex- 
posure, in particular in the outer Galaxy, will be required 
to better determine the longitude profile. 



Using the first year of SPI observations, the emission 
from the central region of the Galaxy could be equally 
well described by various astrophysical stellar bulge and 
halo distributions, or simple analytical shapes [20]. With 
the current, increased data set this is no longer true. The 
preferred stellar bulge models derived by Dwek et al. 
||8[] and Freudenreich [9] from COBE observations in the 
infra-red provide a poorer fit of the 511 keV line distribu- 
tion than simple analytical functions. It therefore appears 
that the distribution of the 511 keV line emission does 
not follow the stellar bulge. Of all astrophysical distribu- 
tions that we tested, only the stellar halo model by Robin 
et al. B3411 could compete. However, given the lack of an 
independent detection of genuine halo emission from out- 
side the bulge region and the inhomogeneity of the expo- 
sure, it is possible that the halo model fares well mainly 
because it provides a good description of the brightest 
emission close to the Galactic center. The good quality 
of the halo model fit, and the results from attempts at fit- 
ting the bulge emission with more than a single ellipsoid, 
indicate that within a few degrees from the Galactic cen- 
ter the annihilation emission exhibits some sort of "peak" 
or "cusp". More observations will be needed to further 
investigate whether the emerging sub-structure of anni- 
hilation radiation from the complex bulge region can be 
associated with gas and/or source distributions. The flux 
attributed to the stellar halo is consistent with the earlier 
analysis of Knodlseder et al. B20I1 and also with estimates 
of halo fluxes by OSSE/SMM/TGRS IMIIl- 

The analyses reported here are ongoing; more refined re- 
sults will be presented elsewhere, including searches for 
point sources. Looking at the years ahead, reducing the 
inhomogeneity of the exposure of the central region of 
the Galaxy is crucial for studying the large-scale distribu- 
tion of annihilation radiation. More observations around 
the Galactic bulge region in particular, but also above and 
below the disk, will be needed to alleviate the currently 
existing systematic uncertainties. 
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Model Component Flux [ph cm 2 s 1 ] 



2 Nested Shells 



0-0.5 kpc Shell 
0.5-1.5 kpc Shell 



(3.88 ±0 
(5.72 ±0 



34) x 10~ 4 
47) x 10~ 4 



3 Nested Shells 



0-0.5 kpc Shell 
0.5-1.5 kpc Shell 
1.5-5 kpc Shell 



(3.88 ±0 
(4.64 ± 
(8.05 ±1 



.34) x 10~ 4 
.52) x 10~ 4 
.76) x IP" 4 



4 Nested Shells 



0-0.5 kpc Shell 
0.5-1.5 kpc Shell 
1.5-5 kpc Shell 
5-8 kpc Shell 



(3.88 ± 
(4.92 ± 
(4.73 ± 1 
(13.6 ±5 



.34) x 10" 4 
.53) x 10~ 4 
.23) x 10~ 4 
.04) x 10~ 4 



2 Nested Shells and Young Disk 



0-0.5 kpc Shell 
0.5-1.5 kpc Shell 
Young Disk 



(3.62 ±0 
(4.74 ± 
(1.11 ±0 



.34) x 10" 4 
.49) x 10~ 4 
18) x 10~ 3 



3 Nested Shells and Young Disk 



0-0.5 kpc Shell 
0.5-1.5 kpc Shell 
1.5-5 kpc Shell 
Young Disk 



(3.66 ± 
(4.21 ±0 
(5.09 ± 1 
(0.94 ±0 



.34) x 10" 4 
.53) x 10~ 4 
.86) x 10~ 4 
19) x 10~ 3 



4 Nested Shells and Young Disk 



0-0.5 kpc Shell 
0.5-1.5 kpc Shell 
1.5-5 kpc Shell 
5-8 kpc Shell 
Young Disk 



(3.68 ±0 
(4.36 ± 
(4.86 ± 1 
(5.82 ± 5 
(0.86 ±0 



.34) x 10~ 4 
.55) x 10~ 4 
.87) x 10~ 4 
.34) x 10~ 4 
.21) x 10~ 3 



2 Nested Shells and Old Disk 



0-0.5 kpc Shell 
0.5-1.5 kpc Shell 
Old Disk 



(3.72 ±0. 
(4.33 ±0. 
(1.82 ±0. 



34) x 
52) x 
29) x 



HP 3 " 

10- 4 

10" 4 



3 Nested Shells and Old Disk 



0-0.5 kpc Shell 
0.5-1.5 kpc Shell 
1.5-5 kpc Shell 
Old Disk 



(3.76 ±0. 
(3.97 ±0. 
(4.33 ±1. 
(1.52 ±0. 



34) x 
54) x 
93) x 
32) x 



HP 3- 

io- 4 
io- 4 

IO" 3 



4 Nested Shells and Old Disk 



0-0.5 kpc Shell 
0.5-1.5 kpc Shell 
1.5-5 kpc Shell 
5-8 kpc Shell 
Old Disk 



(3.75 ±0. 
(4.11 ±0. 
(4.26 ±1. 
(4.35 ± 5. 
(1.41 ±0. 



34) x 
57) x 
93) x 
55) x 

35) x 



10- 4 
io- 4 

10~ 4 
IO" 3 



A 

1273.8 



1294.7 



1302.1 



1310.7 



1318.2 



1319.3 



1312.4 



1317.4 



1318.0 



Table 1. The results of fits in which bulge and/or halo 
emission was modelled by nested shells of homogeneous 
emissivity, and disk emission was described by either the 
young or the old stellar disk model of Robin et al. $3^1 . 
A is the maximum likelihood ratio of the fits. 



